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Abstract 

We study the variational problem that arises from consideration of large deviations for semimartingale 
reflected Brownian motion (SRBM) in R+. Due to the difficulty of the general problem, we consider the 
case in which the SRBM has rotationally symmetric parameters. In this case, we are able to obtain 
conditions under which the optimal solutions to the variational problem are paths that are gradual 
(moving through faces of strictly increasing dimension) or that spiral around the boundary of the octant. 
Furthermore, these results allow us to provide an example for which it can verified that a spiral path is 
optimal. For rotationally symmetric SRBM's, our results facilitate the simplification of computational 
methods for determining optimal solutions to variational problems and give insight into large deviations 
behavior of these processes. 



1 Introduction and Main Results 

In this paper, we analyze the variational problem associated with the large deviations principle for semi- 
martingale reflected Brownian motion (SRBM) in the octant. The SRBM processes of interest arise from 
heavy traffic limits of queueing network processes. Understanding the tail asymptotics of the SRBM's can 
aid in computing their stationary distribution, which in turn gives insight into the behavior of the pre-limit 
queueing processes. 

The typical analysis of large deviations for any process can often be divided into two steps: (1) proving 
a large deviations principle (LDP) and (2) analyzing the resulting variational problem. For particularly 
complex variational problems, one might further subdivide step (2) into: (2a) characterizing optimal paths 
and (2b) optimal path computations. Our primary interest in this paper is in step (2a), especially for SRBM's 
m R%. To gain understanding of the difficulties of the overall investigation, we briefly review some previous 
results. First, with respect to step (1), an LDP for SRBM's in Rl has only been established for special 
cases. For a general dimension d, Majewski examined the special cases of SRBM's arising from feed-forward 
queueing networks [T3] and SRBM's whose reflection matrix is an M-matrix (the so-called Harrison-Reiman 
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case) [TJj. Dupuis and Ramanan [7] obtained an LDP for a generalization of the Harrison-Reiman case. 
It should be noted that these results still leave the LDP for d = 2 unresolved for some parameter cases 
(see for a summary). More recently, Dai and Miyazawa [5] obtained exact asymptotics for SRBM in 
two dimensions using moment generating functions and techniques from complex analysis. However, the 
results are limited to asymptotic behavior along a ray of the quadrant. In a related follow-up paper, Dai 
and Miyazawa 6] provide new insights into the results of Avram et al. |il{ and derive exact asymptotics for 
the boundary measures of SRBM in two dimensions. 

The tasks outlined for step (2) are best explained by examining the case in two dimensions. In this setting, 
Avram et al. [1] and Harrison and Hasenbein (9] gave a complete analytical solution to the variational problem 
for any SRBM of interest (e.g., those possessing a stationary distribution). The analysis was carried out 
in a few steps. First, three general properties of optimal paths were established: convexity, scaling, and 
merging. Second, these properties were used to conclude that only three types of optimal paths are possible. 
Finally, these path properties allow the development of a complete algebraic description of the optimal paths 
in two dimensions. Unfortunately, the situation in three dimensions is considerably more difficult. While 
the properties of convexity, scaling, and merging still apply, they are nowhere near sufficient to characterize 
the possible optimal paths. In order to attack the higher dimensional problem, we examine a special set of 
SRBM cases and develop new techniques for restricting the types of paths which must be examined. 

More specifically, in this paper we investigate the variational problem associated with SRBM in the posi- 
tive orthant for d = 3 in the case in which the SRBM has either rotationally symmetric or mirror symmetric 
data (the latter is a special case of the former) . Note that neither symmetry case we analyze coincides with 
the much studied case of skew- symmetric SRBM's, which have tractable product form stationary distribu- 
tions. 

Our first contribution is to use the Bramson, Dai, and Harrison [3] stability results to derive an appealingly 
simple set of stability conditions for rotationally symmetric SRBM (see Theorem [9] in Section [5]). However, 
the main contribution of the paper is to clarify the nature of optimal paths in three-dimensional variational 
problems and to provide new techniques to achieve this analysis. To best elucidate our contribution, we 
present our main result now: 

Theorem 1. Consider a rotationally symmetric variational problem, as given in Definition^ arising from 
SRBM in R^_. Suppose T = I and 9 < 0. Under Conditions 1 and 2 in Section [$| there always exists 
an optimal path which is either (a) a gradual path (a path which moves through faces of strictly increasing 
dimension) or (b) a classic spiral path. 



Theorem |I9| establishes part (a) and Theorem 21 establishes part (b). An important consequence of this 
result is the following: 

Corollary 2. For the variational problem arising from SRBM in the octant, there exists an example of an 
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optimal spiral path. 

This follows from Theorem [l] and the calculations in Section 10 To the best of our knowledge, this is the 
first time a spiral path has been shown to be optimal for this type of variational problem. Conditions 1 and 2 
are somewhat complicated and we discuss them in detail later. The most important restriction they impose 
is that the SRBM must have reflection vectors which point "outward." We believe that this condition, along 
with the covariance matrix condition, can be relaxed. Note that the negativity condition on the drift 6 is 
not restrictive, since it is equivalent in our case to requiring stability of the associated SRBM. 

An implication of these results is that they provide the basis for tractable numerical methods for comput- 
ing optimal paths in three dimensions. Complementary to our work is a recent paper by El Kharroubi et al. 
[8], which provides some algebraic results for paths in three dimensions. However, most of the results in [8] 
require a priori elimination of certain optimal path types, which we are able to provide. Important related 
computational methodology appeared in Majewski [14] . If one fixes the maximum number of segments in 
the search for an optimal path, Majewski's branch-and-bound algorithm can efficiently produce the desired 
path. 

Our hope is that the path properties we establish can be extended beyond the symmetry cases. However, 
it should be noted that the general case in three dimensions is already known to be fairly complex and in 
fact our main results do not hold for all parameter cases in three dimensions. In [7], the authors show that 
an optimal path to a point in the interior of the octant may have up to five linear pieces, implying that a 
simple characterization of paths in the general d — 3 case is non-trivial. This five-piece path is depicted in 
Figure [l] (in the fi gures in this paper, dotted lines indicate a segment contained in the interior of the octant). 
Nonetheless, we believe that deriving new properties of optimal paths for special cases is a necessary building 
block for solving other variational problems, and provides for a better understanding of the tail asymptotics 
for SRBM. 

This paper is structured as follows. In Section [2] we introduce the variational problem (VP) analyzed 
throughout the paper. This problem arises from studying large deviations of SRBM in the orthant, concepts 
that are described in Section [3j Section [4] introduces the symmetric cases of the SRBM and VP which are 
of interest to us. In Section [5j we use the framework of Bramson et al. [3] to derive the stability conditions 
of symmetric SRBM. The next three sections characterize the nature of (piecewise linear) optimal paths 
with a finite number of segments. Section [9] is devoted to paths with an infinite number of pieces and it is 
demonstrated that only classic spiral paths can be optimal. Finally, in Section [10] we provide an example of 
a spiral path that is indeed optimal. 
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Figure 1: An Optimal Path with Five Pieces 

2 The Variational Problem 

In this section, we define the variational problem of interest in this paper. First, we give notation and 
definitions which follow as closely as possible to those given in Avram et al. PQ. 

Let d > 1 be an integer and 9 a constant vector in R d . Also, r is a d x d symmetric and strictly positive 
definite matrix, and R is a d x d matrix. The triple (9, T, R) provides the data to variational problems 
and, as described later, associated reflected Brownian motion processes. Throughout the paper, all vector 
inequalities should be interpreted componentwise and all vectors are assumed to be column vectors. Finally, 
for vectors v € R d and w € R d we define the inner product 

(v, w) = v'T~ 1 w 

and the associated norm ||u|| = (v, v). 

In order to more easily define the VP, we first introduce the Skorohod problem associated with the 
matrix R. Thus, let C([0, oo),H d ) be the set of continuous functions x : t G [0, oo) — > x(t) € R d . A function 
x G C([0, oo), R d ) is called a path and is often denoted by x(-). We now define the Skorohod problem 
associated with a reflection matrix R. 

Definition 1 (The Skorohod Problem). Let x be a path. An R-regulation of x is a pair of paths (z,y) £ 
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C([0,oo),R d ) x C([0,oo),R d ) such that 



z{t) 



x(t) + Ry(t), t>0, 



(1) 



z{t) > 0, t > 0, 



(2) 



y(-) is non-decreasing, j/(0) 



(3) 



/ Zi(s)dyi(s) = 0, i=l, 



,rf. 



(4) 



Jo 



When the i?-regulation (y, z) of a; is unique for each x £ C([0, oo),R ), the mapping 



-0 : a; — > 0(a;) 



is called the reflection mapping from C([0,oo),R d ) to C([0, oo), R+ ). When the triple (x.j/.z) is used, it is 
implicitly assumed that (y, z) is an i?-regulation of x. 

An important issue when defining the Skorohod problem is whether a solution exists for any given path 
x. If the reflection matrix R is completely-5, as defined below, then indeed there is a solution for every x 
with x(0) > (see Bernard and El Kharroubi [2]). 

Definition 2. A d x d matrix R is said to be an 5-matrix if there exists a u > such that Ru > 0. The 
matrix R is completely-5 if each principal submatrix of R is an S-matrix. 

The class of "P-matriccs, defined below, also plays an important role in the development of SRBM theory 
and associated variational problems. 

Definition 3. A d x d matrix R is said to be a "P-matrix if all of its principal minors are positive. 

In addition to the issue of existence of solutions to the Skorohod problem, there is also the matter of 
the uniqueness of the solution, for a given path x. It is useful when defining the VP to have a notational 
convention which applies when solutions are not unique. To this end, we assume that if the Skorohod 
problem is non- unique, then tp{x) represents a set of paths (solutions) corresponding to x. Furthermore, the 
expression 



1>(x)(T) 



indicates that there exists az£ ip{x) such that z(T) = v. 

We now define the variational problem studied in this paper. 



Definition 4 (The Variational Problem). 




(5) 



where % d is the space of all absolutely continuous functions x(-) : [0, oo) — >■ R d which have square integrable 
derivatives on bounded intervals and have x(0) = 0. 
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Definition 5. Let v G R + . If a given triple of paths (x,y,z) is such that the triple satisfies the Skorohod 
problem, z(T) = v for some T > 0, and 



then we will call (x, y, z) an optimal triple, for VP |5|), with optimal value I(v). The function z is called an 
optimal path if it is the last member of an optimal triple. Such a triple (x, y, z) is also sometimes referred 
to as a solution to the VP (^J). T is called the optimal time for such a solution. 

3 SRBM and Large Deviations Background 
3.1 Semi-martingale Reflected Brownian Motion 

We now define the semi-martingale reflected Brownian motion (SRBM) on the positive orthant associated 
with the data (9,T,R). Let B denotes the c-algebra of Borel subsets of Hi.. A triple (Q,J 7 ,{J r t}) is called 
a filtered space if il is a set, J 7 is a cr-field of subsets of f2, and {J 7 *} = {.Ft, t > 0} is an increasing family of 
sub-cr-fields of J 7 , i.e., a filtration. If, in addition, P is a probability measure on (f2, J 7 ), then (0, J 7 , {J~t}, P) 
is called a filtered probability space. 

Definition 6 (SRBM). Given a probability measure v on (R^,£>), a semi-martingale reflecting Brownian 
motion associated with the data (9,T, R,v) is an {J^t}- adapted, d-dimensional process Z defined on some 
filtered probability space (Q, {jF t },Pv) such that 

(i) P„-a.s., Z has continuous paths and Z(t) € R+ for all t > 0, 



(a) X is a d-dimensional Brownian motion with drift vector 9, covariance matrix Y and X(Q) has 
distribution v , 




(ii) Z = X + RY, P v -a.s., 



(Hi) under P 



(b) {X(t) - X(0) -9t,F t ,t> 0} is a martingale, 



(iv) Y is an {Tt}- adapted, d-dimensional process such that P„-a.s. for each j = 1, . . . , d, 



(a) Y J (0)=0, 




G 



An SRBM associated with the data (R/j_, 9, T, R) is an {Tt} -adapted, d- dimensional process Z together with 
a family of probability measures {P X7 xeR d + } defined on some filtered space (CI, J- , {J 7 *}) such that, for each 
x G R+, (i)-(iv) hold with P u = P x and v being the point distribution at x. 

Recall that the parameters 9, V and R are called the drift vector, covariance matrix and reflection matrix 
of the SRBM, respectively. The results of Reiman and Williams [15] and Taylor and Williams [16] imply 
that the necessary and sufficient conditions for the existence the SRBM is that R is completely-5. 

The measure v on (R+, B) is a stationary distribution for an SRBM Z if for each A G B, 

v(A) = j P x {Z(t) G A} v{dx) for each t > 0. (6) 

When v is a stationary distribution, the process Z is stationary under the probability measure P„. In our 
discussion below, we are concerned only with the (unique) stationary distribution for the SRBM with data 
(9, R, r) and therefore we drop the v notation. 

3.2 Large Deviations 

The motivation for studying the variational problem introduced in Section [2] comes from the theory of large 
deviations. For SRBM's in R^, we have the following statement of the large deviations principle, which has 
only been established for some special cases, as noted in the introduction. 

Conjecture 3 (General Large Deviations Principle). Consider an SRBM Z with data (9,T,R). Suppose 
that R is a completely-S matrix and that there exists a probability measure under which Z is stationary. 
Then for every measurable A C R^ 

limsup — logP 7r (Z(0)/M G A) < - inf I(v) (7) 
u->oo u veA? 

and 

liminf - logP^(Z(0)/u G A) > - inf I(v) (8) 

u— >oo u VGA" 

where A c and A° are respectively the closure and interior of A. 

The specific connection between this LDP statement and the VP is that the function /(•) appearing 
above is the same function which appears in Definition [5] 

4 Symmetric SRBM 

In this paper, we study solutions to the variational problem associated with the LDP introduced in the 
previous section. The three-dimensional case is considerably more difficult than the two-dimensional case 
and thus we confine our study to SRBM's with some symmetry in the data. The special cases we study 
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are called rotationally symmetric and mirror symmetric and are denned below. These symmetries, while 
restrictive, provide a considerable simplification of the analysis. 



Definition 7. For d = 3 the data (9,T,R) is 
three conditions hold: 

1. R has the form 



R = 



to be rotationally symmetric if all of the following 



I I r 2 r± \ 

n 1 r 2 
\r 2 n 1 J 



2. The drift has the form 9 = (9q, 9q, 9q)' . 

3. The covariance matrix has the form 



r 



f o 2 pa 2 pa 2 



pa a 



y pa 2 pa 2 a 2 ) 



where — 1 < p < 1. 



Some statements in the rest of the paper relate only to R and in this case we call R alone rotationally 
symmetric if and only if R has the form given in the definition above. We employ a similar convention for T. 

Definition 8. For d = 3 the data (9, T, R) is said to be mirror symmetric if it is rotationally symmetric 
and in addition r\ =r 2 . 

Notice that the rotationally symmetric T matrix also appears to be mirror symmetric. Since covariance 
matrices are by definition symmetric (in the standard matrix algebra sense), there is no sensible way to define 
a rotationally symmetric T which is not also mirror symmetric. For a rotationally symmetric T we have the 
following result, proved in the Appendix, which will be used in demonstrating optimal path properties. 

Lemma 4. If T is rotationally symmetric, then 

' 7o 7i 7i ^ 
1 -a 71 70 7i 

\ 7i 7i 7o j 
with 70 > 7i- 

Some readers may also be familiar with the skew-symmetry condition (see [101 111] ) which is 

2T = RD~ 1 A + AD^R', (9) 
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Is R completely-S? 



SRBM not defined 



isi? _1 e < o ? 



SRBM not stable 



IS (6»,R) 6 C ? 



Is there a divergent 
solution of the LCP? 



is p(e,R) < 1? 



SRBM not stable 



SRBM stable 



SRBM stable 



Figure 2: Existence and Stability of SRBM in the Octant 

where D — diag(i?) and A = diag(r). This condition is necessary and sufficient for the stationary density 
function of the SRBM to admit a separable, exponential form. Our notions of symmetry do not coincide 
in any meaningful way with the notion of skew-symmetry. It can be checked that rotationally symmetric 
SRBM data is also skew-symmetric if and only if r\ + r2 = 2p. 

In subsequent sections, we provide results for both SRBM and the associated variational problems. Thus, 
for an SRBM with rotationally symmetric data we use the abbreviation RS-SRBM. Similarly, for an SRBM 
with mirror symmetric data we use MS-SRBM. The associated variational problems take the same data and 
when stating results for VPs we use the abbreviations RSVP and MSVP. 

5 SRBM Stability Conditions 

For SRBM in three dimensions Bramson et al. [3] obtained results which, in addition to previous results, 
give a complete characterization of existence and stability of SRBM. This characterization is summarized in 
Figure [2J The results of this section specialize their results for RS-SRBM. First, however, we need to define 
a few terms appearing in the figure. 

We define the solutions to the linear complementarity problem (LCP) in dimension d as follows (see |3] 
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for background). Vectors u,v <E R d comprise a solution to the LCP if 



u, v > 

v = 9 + Ru 
u ■ v = 0. 

Using the terminology in [3J a solution (u, v) to the LCP is called stable if v — and the solution is called 
divergent otherwise. The existence or non-existence of a solution to the LCP must be checked in the bottom 
decision point in Figure [2] 

Bramson et al. [3J also define various subsets of the data pairs (9, R) which relate to the third line 
of decision points in Figure [2] To avoid overlapping notation we specialize their definitions now to the 
RS-SRBM case. 



First, for a pair (9,R), 



C\ = {(0 o ,ri,r 2 ) : 9 < 0,r x < l,r 2 > 1} and 
C 2 = {(0o,ri,r 2 ) : 9 < 0,n > l,r 2 < 1}, 



with C = C 1 \JC 2 . 

Next, for a data pair (9, R) € C\, 

and for (9,R)e C 2 



,R) 



13(9, R) 



, 3 

r 2 



n-l 



1 \ 3 

ri - 1 



For general SRBM data /3(0, i?) depends on 9 but in the rotationally symmetric case the dependence disap- 
pears. These definitions are related to spiral piecewise linear solutions of the Skorohod problem. Section 3 
in [3j should be consulted for an in-depth explanation of how these expressions arise. 

We are now prepared to present a series of lemmas which lead to the main stability result of this section. 
The first lemma probably appears in a textbook somewhere, but we state it here and prove it in the Appendix 
for completeness. For later use, note that the lemma implies that a + 6 + c ^ 0. All the results stated in this 
section apply to the three-dimensional case. 

Lemma 5. If a reflection matrix R is non-singular and rotationally symmetric then its inverse must be of 
the form 

1 a b c\ 
R 1 = cab 
^ b c a 

and (a + b + c)(l + r x + r 2 ) = 1. 
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Since existence of an SRBM requires that R is completely-^, our first task is derive a simple condition 
to insure that this characterization holds. 



Lemma 6. Suppose a matrix R is rotationally symmetric. Then R being completely-S is equivalent to 
1 + r x + r 2 > 0. 

Proof. First suppose R is completely-5 and l + rx+r 2 < 0. We derive a contradiction. Since R is completely- 
S, there exists a vector u = (ui, u 2 , u^)' > such that Ru > 0. Summing the equations in Ru > we have 



(1 + r\ + r 2 )(u 1 + u 2 + u 3 ) > 0. 



(10) 



But if 1 + ri + r 2 < 0, then there is no u > satisfying (10 1, which is a contradiction. So, we have proven 



that completely-5 implies 1 + n + r 2 > which is one direction of the equivalence. 

Now assume that 1 + r\ + r 2 > 0. Then note that u = (1, 1, 1)' satisfies Ru > 0. This implies that R 
is an 5-matrix. We must now verify that the two-by-two principal submatrices are also 5-matrices. These 
submatrices take the form 



Si = 



1 r 2 

n l 



and S 2 = 



i n 

ri 1 



We prove the result for S\ only, since the argument for 5*2 is completely analogous. Now, for S\ to be an 
5-matrix there must exist an (1*1,1*2)' > such that, 



U\ + u 2 r 2 > 

U\T\ + U 2 > 0. 



(11) 

(12) 



First suppose ri,r 2 > 0. Then any (ui,u 2 )' > satisfies (11) and (12). In the cases (i) r% > 0, r 2 < 0, 



(ii) ri < 0, r 2 > 0, and (iii) ri,r 2 < 0, it can be checked that (11) and (12) holding for any positive it is 
equivalent to r\r 2 < 1. This last inequality obviously holds in cases (i) and (ii). It also holds in case (iii) 
because 1 + r\ + r% > insures r l5 r 2 > —1 and at least one of them is strictly bigger than -1. Finally, if 



ri = and/or r 2 = then, for example u — (1, 1) satisfies (11) and (12) 



□ 

Having now dispatched with the first line in Figure [2] we present a lemma relating to the second line. 

Lemma 7. Let the data (9,T,R) be rotationally symmetric and let R be non-singular and completely-S. 
Then R _1 6 < is equivalent to 9 < 0. 

Proof. Using Lemma [5] the condition R~ 1 9 < reduces to 

(a + b + c)9 Q < 0. 
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The second part of Lemma [5] states that (a + b + c) = [1 + r% + r 2 ] 1 - Hence, we can rewrite the condition as 

R- X e = < 0. (13) 

1 + n + r 2 

By Lemma [6j the completely-iS condition is equivalent to 1 + r\ + r 2 > 0. Given this inequality, (131 is 

clearly equivalent to 9q < 0. □ 

Now we proceed to results involving the last two lines in Figure [2j 

Lemma 8. Let the data (9, T, R) be rotationally symmetric. Suppose further that R is non-singular, 
completely-S , and R~ 1 < 0. Then the SRBM associated with (9,T,R) is stable iff ri + r 2 < 2. 

Proof. Our proof relies on the results in [3] as depicted in Figure [2] The assumptions of the lemma place us 
in the lower half of the figure. To further partition the proof, we divide the (ri,?^) plane into four regions: 

• C x ={(ri,r 2 ) eH 2 :n < l,r 2 > 1} 

• C 2 = {(r 1 ,r 2 )eH 2 :r 1 >l,r 2 <l} 

• C 3 = {(n, r 2 ) G H 2 : r x > 1, r 2 > 1} \ (1, 1) 

• C A = {(n,r 2 ) G H 2 : n < 1, r 2 < 1} \ (1, 1). 

We do not include the completely-5 condition that r\ + r 2 > —1 in this partitioning scheme because the 
condition is not employed directly in the algebraic arguments below. Under our assumption R~ x 9 < 0, the 
definitions of C\ and C 2 coincide with the Bramson et al. [5] definitions given in Section [5j Furthermore, 
note that if (n, r 2 ) G C x U Ca then (0, i?) G C. 

The one point of the plane not included in the union of these regions is r*i = r 2 = 1. The matrix R is 
singular in this case, which violates the assumption of the lemma. Notice that the line r\ + r 2 = 2 bisects 
C\ 1J C 2 and that C3 lies entirely above this line and C4 entirely below this line. 

Case 1: Suppose (ri,r 2 ) G CilJC 2 . In this case, stability of the SRBM is equivalent to (3(9, R) < 1. Now, 
when (ri,r 2 ) G C± we have, 

P(9,R)=' ' " " 



Since the numerator and denominator are both negative for (r l5 r 2 ) G C\, the condition fi(9,R) < 1 is 
equivalent to 1 — r 2 > r% — 1, which holds iff ri + r 2 < 2. 
Next, when (ri,r 2 ) G C 2 , 

/3(0,i?) = 



1 - n 



J2-1, 

Again, the numerator and denominator in the last expression are both negative for (ri,r 2 ) G C 2 . Therefore, 
(3(9, R) < 1 is equivalent 1 — r\ > r 2 — 1, which also holds iff ri + r 2 < 2. 
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So, for all of Case 1, r± + r 2 < 2 is necessary and sufficient for stability. 

Case 2: Suppose (n, r 2 ) € C3 [J C4. In this case stability of the SRBM is equivalent to the nonexistence of 
a divergent solution to the LCP. 

Case 2a: We examine the case (ri, r%) <E C3. Since ri + r 2 > 2 for all points in C3 we need to show instability 
for data in this region. In the LCP, let u = (—60, 0, 0)' and v = (0, r\ — 1, r 2 — 1)'. This is clearly a divergent 
solution to the LCP for any (ri,r2) € C3. Therefore the corresponding SRBM is never stable in this case. 

Case 2b: We examine the case (r l7 r 2 ) € C4. Since r x + r 2 < 2 for all points in C3 we need to show stability 
for data in this region. 

Consider a solution u, v to the LCP. We show that there exist no divergent solutions for this case. If u > 0, 
then v = and the solution is stable. If u = 0, then we must have v = 9 < which is not an allowable LCP 
solution. So if there exists a divergent solution, either one term or two terms in u = (iti, u 2 , U3)' is positive. 
Suppose one term is positive and it is u\, which implies v\ = 0. Then we have v\ = 9q + U\ = yielding 
iti = —9q. Therefore, in this case the unique a solution to the LCP must be of the form u = (— 60, 0,0)', 
v = (0, ?'i — 1, r- 2 — 1)', which violates the non-negativity condition of v. Exactly analogous arguments show 
that neither u 2 nor 113 can be the positive term. So, there exist no LCP solutions in which only one term in 
u is positive. 

Next, suppose that two terms of u are positive. Again, without loss of generality, suppose u\ > 0, w 2 > 0, 
and W3 = 0, which implies v\ = w 2 = 0. From the LCP equations we have u\ + r 2 w 2 + 9q — and 
u 2 + riui + 9 = 0. Solving these yields: 

1 - rir 2 1 - rir 2 

If ri = 1, then r 2 < 1 and these equations force u 2 = 0, which contradicts our assumption on u. Similarly, 
we cannot have r 2 = 1. So, we now assume that both r\ and r 2 are strictly less than 1. In this case, the 



solution given in ( 14 1 implies that both u\ and u 2 are positive. Once again using the LCP equations we 
obtain: 

ri - r\ + r 2 - r\ + r x r 2 - 1 
«3 = -fo- : ■ (15) 

1 - nr 2 

Note that 1 - r x r 2 > 0, -6 > 0, and 

r\ - r\ + r 2 - r\ + rir 2 - 1 < r\ + r<i — r x r 2 - 1 = -(1 - ri)(l - r 2 ) < 0, 

Therefore, W3 < which violates the non- negativity condition in the LCP. 

We have now demonstrated that no divergent LCP solutions exist when (r\,r 2 ) S C4. So, any SRBM 
with data in this region is stable. □ 

Lemmas [6] through [8] then imply simple existence and stability conditions for RS-SRBM in three dimen- 
sions. 
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Theorem 9. Consider an SRBM in three dimensions with rotationally symmetric data (9,T,R). The 
necessary and sufficient conditions for existence and stability of such an SRBM are 9q < and — 1 < 
n + r 2 < 2. 

The results in [5], which we shall make use of in later sections, require that R be a "P-matrix. The next 
result shows that this is not a restriction in the rotationally symmetric case, given that we only study stable 
SRBM's. The proof is given in the Appendix. 

Theorem 10. Suppose R is rotationally symmetric and ri+r 2 < 2. Then R being completely-S is equivalent 
to R being a V -matrix. 

6 Optimal Path Preliminaries 

In this section we establish some notation and review the properties pertaining to optimal paths. As much 
as possible we use notation which is consistent with either [1] or [8] . Many of our results rely on algebraic 
expressions given in [3]. First, we give expressions for the optimal costs of various types of paths. 
Set I = {1, 2, . . . , d} and for K C I define the face associated with K as follows: 



F K = {v G : Vi = for all i € K). 



When d — 3, if \K\ = 2 then Fk is an axis and if \K\ = 1 then F% is a 2-dimensional face. 



Definition 9. We define the following costs. 



1. (Direct Path Cost) For »el|, set 




2. (One-piece Reflected Path Cost) For K <Z I and v € Fk, set 




3. (Two- Piece Path via Face F K ) Let d = 3 and \K\ < 2. For v € set 




inf (X K (w) +l (v - w)). 



4-. (Two-Piece Path via an Axis). Let d = 3 and \K\ = 2. For i £ K and v € Fi, set 



^K,i ( v ) = inf (Ik(w)+1 1 (v-w)). 
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5. (Three-Piece Gradual Escape Path) Let d = 3, \K\ = 2. For i E K and v £ R+, set 

Z*»= inf (X 2 Kti (u)+Xo (v-u)). 

Each cost above corresponds to the cost for an optimal path of a certain type, as denoted in each item 
in the list. These costs, and the associated paths, are the building blocks for constructing paths which are 
optimal in the original variational problem. 

In PP, the authors established various properties of optimal paths that hold in all dimensions. The first 



three items in Lemma 11 restate those properties. We add a fourth property for RSVPs and MS VPs in three 
dimensions. These properties are frequently used to establish results in subsequent sections. The first three 
parts are proved in pQ, the fourth result is evident using symmetry. In the Appendix, we state and prove 



a simple extension to the convexity property (see Lemma 23 1 . The convexity property below implies that 
direct paths within a face should have constant velocity and direction. The extension shows that this also is 
true for reflected paths. 

Lemma 11. 1. (Optimality of Linear Paths via Convexity) Let g be a convex function onM. d , andx € T-L d . 
Then for t\ < ti, 

This implies that a (one-piece) linear path minimizes the unconstrained variational problem. 

2. (Scaling) Consider a variational problem with v 6 Mi. ForVk > 0, L(kv) = kl(v). Furthermore, 
if (x,y,z) is the optimal triple for v and x,y,z is the optimal triple for kv, then x(t) = kx{t/k), 
y(t) = ky(t/k), z(t) = kz(t/k). 

3. (Merging Paths) Let {x\,yx,zx) be an R-regulation triple on [0,ti] with Zi(0) = and z\{ti) = w and 
(x2,2/2, 3a) ^ e an optimal triple on [32,^2] with 22(52) = w and 22^2) == v - Suppose both X\ and X2 are 
absolutely continuous. Define 



z{t) = 



x(t) 



«l(t) forO<t<t u 

z 2 (t — + s 2 ) for t\ < t < ti + ti — s 2 , 

x 1 (t) forO<t<tt, 

X2(t — ti + s 2 ) for ti < t < tx + t 2 - S2, 



y(t) - 



Vi{t) forO<t<t u 
Viit -ti + s 2 ) for t\ <t <h+t 2 - s 2 , 
and s = t\ + t% — S2- Then (x, y, z) is an R-regulation triple on [0, s] with z(0) — and z(s) = v. 



15 



4- (Symmetric Terminal Points) Consider an RSVP. If v± = (a,b,c), v 2 = (6, c, a) 7 V3 — (c, a, 6) with 
a, 6, c > 0, then I(v\) = I{v2) = -^(1*3). Furthermore, each optimal path to one of these points is 
rotationally symmetric translation of an optimal path to one of the other points. For an MSVP case, 
let v[ = (6, a, c), v 2 = (c, b, a), v' 3 = (a, c, b), then I(vi) = I(v 2 ) = I(v 3 ) = I(v[) = I(v' 2 ) = I(v' 3 ). 



7 Eliminating Bad Faces in RSVPs 

The next result is one of the key results in the paper, since it allows us to eliminate entire categories of paths 
by eliminating paths whose penultimate pivot point is on a "bad face." Below, the distance between a face 
and a point is the standard Euclidean distance from a point to the associated face. 

Theorem 12 (Bad Faces). Consider a variational problem with terminal point v € int(R + ) and consider 
an optimal triple (x, y, z) to v. Let w be the last point of z which is not in int(R 3 ^). Then 

(a) If the variational problem is an RSVP, then there exists an optimal path for which w is in one of the 
two nearest faces to v. 

(b) If the variational problem is an MSVP, then there exists an optimal path for which w is in the nearest 
face to v. 

Proof. Let the terminal point be v = (v±, v 2 , V3), and without loss of generality, assume v 3 > v 2 > v± > 0. 
Define u 1 = (0, a, b), u 2 = (b,0,a), u 3 = (a, 6,0) and w 4 = (0,6, a). Furthermore, we assume that we cannot 
have both a — and 6 = 0. 

First, we want to compare the optimal cost from u l to v for various values of i. Of course, by convexity 



(Lemma 11 part 1), the optimal path from any u l to v must be a linear path. Now, recall that 

2b(«-u*) = \\e\\\\v - u l \\ - (9,v-u l ), 

for i G {1, 2, 3, 4}. It can be checked that 

(9, v - u l ) = 6 a cr~ 2 (2ji + 7o)(vi +V2 + V3- a-b). 

Hence this portion of the cost is independent of i. So it is sufficient to analyze \\v — u l \\ or, equivalently, 
\\v-u l \\ 2 . Note then that 

||v - u l \\ 2 = {v,v) + (u\ u l ) - 2{v, u l ), 

where (u l ,u l ) = a~ 2 [(a 2 + b 2 )-fo + 20671]. Therefore, the first two terms in \\v — u l \\ 2 are also independent 
of i. So, finally, we have 

(v^u 1 - u 3 ) = er~ 2 [(7o - 7i)(t>3 - v 2 )b + (70 - 7i)(«2 - Vi)a] > 0, 



16 



where the inequality is due to our assumption on v and Lemma|4] This of course implies that (v, u 1 ) > (v, u 3 ) 
with equality iff (V3 — v 2 )b + (v% — Vi)a — 0. Therefore we have 

io{v-u 3 ) yioiv-u 1 ). (16) 

Now, let w be the last point which is not in the interior of the octant, for an optimal path with terminal 
point v. By convexity, the segment wv must be linear. Suppose then that w is in F 3 . If v 3 ^ v 2 then F 3 
is the furthest face from v. Thus w must be of the form of u 3 and accordingly we take w — u 3 — (a, b, 0). 
Now, consider the point u 1 — (0, a, b). Note that the optimal cost from the origin to u 3 and the optimal cost 
from the origin to u 1 must be equal due to rotational symmetry. By the merging and convexity properties 



of Lemma 11 to establish (a) it suffices to show that the optimal cost from u 1 to v using a direct path is 
less than or equal the cost from u 3 to v via a direct path. This result was already demonstrated, as seen in 



(16). 

At this point some discussion may be needed to see that part (a) of the theorem has been proved. Suppose 
first that i>3 > v 2 > v\ > 0. Then F$ is the unique furthest face from v. In this case we can strengthen the 
conclusion of part (a). In particular, the last boundary point in an optimal path must emanate from one 
of the two nearest faces. Next, if v 2 = i>3 then all three faces can be classified as "one of the two nearest" 
and the statement of (a) holds by default. Finally, if v 3 > v 2 = «i > there are two cases. If b 7^ 0, then 
u 3 is in the interior of F3 and there must exist a strictly cheaper path through u 1 . If b = 0, then the cost of 
the paths through u 3 and u 1 are the same. Either path is considered to be via ^(albeit on the boundary) 
which is one of the two nearest faces to v. Hence, the result in (a) is still valid. 

We now address part (b) of the theorem. First, it can be checked that 

{v,u A -u 2 ) = a- 2 [b(v 2 - Vl )( l0 - Ji)} >0, 
with equality if b — or v\ = v 2 . Using the calculations from the RSVP case, we have 

Mv-u 2 ) >i («-« 4 ). (17) 

By mirror symmetry, the optimal cost from the origin to u 3 and the optimal cost from the origin to u 4 must 
be equal. 

The remainder of the proof is similar to the part (a) argument. Again, let w be the last point which is 
not in the interior of the octant, for an optimal path with terminal point v. Suppose first that w is in F 3 . 
Unless vi = v 2 — W3 (in which case the result holds trivially), then F3 is one of the two furthest faces from v. 
Recall that an MSVP is also an RSVP, so we can apply part (a) of the theorem to conclude that there must 
exist an optimal path to v with w € F 2 . Without loss of generality, assume then that w = u 2 = (b, 0, a). By 
mirror symmetry, the optimal cost from the origin to u 3 and the optimal cost from the origin to u 4 must 



be equal. However, by (17) the optimal cost from u to v using a direct path is less than or equal the cost 
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from u 2 to v via a direct path. Hence, there exists a path for which u 4 is last point on the boundary of the 
octant, with lower (or equal) cost to the path through u 3 . 

As in part (a), there are some special cases in part (b), specifically, when V\ = v 2 or b = 0. If V\ = v 2 
then both F\ and F 2 are considered the "nearest face" and the statement holds immediately by applying 
part (a). If v 2 > «i and b =^ 0, then there exists a strictly cheaper path through u 4 . If v 2 > v% and b = 0, 
then u 2 is considered to be in F\ and again the result holds. □ 

The easiest way of rephrasing the RSVP result is as follows. Consider a terminal point v with a unique 
farthest face. Then the last linear segment in an optimal path cannot emanate from the interior of the 
farthest face. Similarly, for an MSVP with a unique nearest face to the terminal point v, the last linear 
segment must emanate from the nearest face. 

Note that the results in Theorem [12] are proved only for v in the interior of the octant. The arguments in 
the proof of the theorem lead immediately to the following extensions for terminal points on the boundary 
of the octant. 



Remark 1. When V\ = and v 2 > 0, v is in the interior of F\. The first part of Theorem 12 holds in the 
following sense: For an RSVP, there exists an optimal path whose last segment does not emanate from the 
interior of F3. Furthermore, the last segment can not originate from F 2 ^ although it may originate from 
F1.3. 



Remark 2. When v\ = v 2 = and V3 > 0, v is on the axis F\ t2 . Again, the first part of Theorem 12 
holds. In particular, there exists an optimal path whose last segment does not emanate from the interior of 
the farthest face, which is F 3 in this case. 

Finally we believe that the results of this section can be generalized to higher dimensional RSVPs and 
MSVPs with minor modifications to the proofs. 



8 Further Optimal Path Characterizations 

Our eventual goal is to show that optimal paths in three dimensions can be of only two types: gradual paths 
and classic spirals. Demonstrating this requires the establishment of a number of properties for paths with 
a finite or infinite number of linear segments. The results in this section are related to paths with a finite 
number of segments, although some of these properties are used later on to establish characterizations for 
paths with an infinite number of segments. 

In various proofs in this section it is useful to consider paths (and the corresponding costs), which are 
feasible, but not necessarily optimal. Therefore, we introduce the following definition. 
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Definition 10. (Cost of a Feasible Path) Given a one-piece feasible R-regulated triple (x,y,z) from u to a 
terminal point v, define the corresponding cost along that path to be 

H y , z {v-u)= l -J*\\x{t)-6\\ 2 dt, 

where z(T) = v — u. 

The next several results provide detailed characterizations of optimal paths. Unfortunately, the overall 
connection will not be apparent until we bring them together to prove the main results. 

Lemma 13 (The Switchback Lemma). Consider a VP with T = I. Let v 1 , v 4 <G int(Fi) and v 2 ,v 3 e int(F 2 ). 
Then the path from v 1 to v 4 consisting of the following linear segments is strictly suboptimal: a direct segment 
from v 1 to v 2 , a reflected segment from v 2 to v 3 , a direct segment from v 3 to v 4 . 

Proof. Let v 1 = (O,^,^), v 2 = (v 2 ,0,v^), v 3 = (v 3 ,0, and v 4 = (0, v 4 ,v 4 ). Assume first that v 2 > v 3 . 
Define is 2 = (v\ — v 3 , 0, v 2 ) and v 3 = (0, 0, v 3 ) which are both in F 2 . Consider a new path from v 1 to v 4 as 
follows: a direct segment from v 1 to v 2 , a reflected segment from v 2 to v 3 , a direct segment from v 3 to v 4 . 
We show that the new path has a strictly lower cost than the original path. Notice that v 3 — v 2 = v 3 — v 2 
so it suffices to compare Io(v 2 — v 1 ) +X (v 4 — v 3 ) with Xo{v 2 — v 1 ) +X (v 4 — v 3 ). By definition 

l {v 2 - v 1 ) +lo(v 4 - v 3 ) = \\9\\(\\v 2 - w 1 !! + ||i; 4 - v 3 \\) - (0,v 2 - v 1 + v 4 - v 3 ) 

and 

l (v 2 - v^+Mv 4 - v 3 ) = \\9\\(\\v 2 -v 1 ^ \\v 4 - v 3 \\) - (6,v 2 - v 1 + v 4 - i 3 ). 
It is easy to check that 

(9, v 2 -v 1 +v 4 - v 3 } = (6, v 2 - v 1 + v 4 - v 3 } 
so it is enough to compare \\v 2 — v 1 ]] + \\v 4 — v 3 \\ with ||u 2 — v 1 ]] + \\v 4 — v 3 \\. Now when r = /, we have 

(\\v 2 -v 1 ^ \\v 4 -v 3 \\)- (\\v 2 - + ||t; 4 -V 3 \\) = ^P+(v 2 ) 2 + ^q+(v 3 r - yfp + (v 2 - v 3 )* - > 

where p = (vi,) 2 + (w| — w|) 2 > and q = (v 4 ) 2 + (v 4 — u|) 2 > 0. Thus, the newly constructed path has 
a strictly lower cost. If v\ < vf then re-define v 2 — (0,0, v 2 ) and v 3 = (vf — vf 7 0,vl). The proof of the 
corresponding result for this case is analogous to the first case. □ 

This result is the key to showing that "exotic" paths which seem intuitively "bad" are indeed suboptimal. 
In particular, it shows that paths which switch back and forth between two faces are not cost effective. 
Analogous arguments show that the lemma holds for any pair of two-dimensional faces. 

The next result is important in establishing the optimality of gradual paths. In this and later proofs, we 
use the standard notation = (0, 0, 1). 
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Lemma 14. Consider an RSVP with V = I and r 2 > 0. Let v 1 — (0,v 2 ,u|) an d v 2 — ( w i>0,i> 3 ), such that 
v 1 G int(Fi) and v 2 G int(F 2 ). Then the path from v 1 to e 3 consisting of the following linear segments is 
strictly suboptimal: a direct segment from v 1 to v 2 and a reflected segment from v 2 to e 3 . 

Proof. Define v 2 = (0,0, vf). We show that 

Mv 2 - v 1 )+l 2 {e 3 - v 2 ) > l (v 2 - v 1 ) +i 2 (e 3 - 5 2 ), 

implying that there exists a better path from v 1 to e 3 , via v 2 . 

Suppose (x 1 jV 1 , z 1 ) is an optimal triple from v 2 to e 3 with corresponding time T 1 , and let (x 2 ,y 2 ,z 2 ) 
be the optimal triple from v 1 to v 2 with corresponding time T 2 (since this path is direct x 2 — z 2 ). Set 
(x 1 ^), ?/(£), z 1 ^)) = (i 1 ,^ 1 ,^ 1 )*, x 2 (t) = x 2 t, with z 1 = (z\,z\,z\)' . Notice that y 1 = (0,^,0)' and 
z\ = 0. Let z 1 (t) = t(0, 0, z\)' and i 1 (t) = t(—r 2 y 2 , — y 2 , x \)' ■ It can be checked that (a; 1 , y 1 , z 1 ) is a feasible 
triple from v 2 to e 3 with T 1 = T 1 . Similarly, setting x 2 (t) = t(0,x 2 ,x 2 )' yields the feasible triple (x 2 ,0,x 2 ) 
from v 1 tov 2 , with f 2 =T 2 . 

Using the paths defined above we have: 

Me 3 - v 2 ) = ^[(-zl - v 2 y\ - d Q ) 2 + {-y\ - O ) 2 + (1 - z\ - r lV \ - 9 Q ) 2 } 



and 



Mes - v 2 ) < hr^i-rtyl - 9 ) 2 + {-y\ - Q ) 2 + (1 - z{ - r lV \ - 6 )% 



where the inequality is due to the fact that (x 1 ,?/ 1 ,^ 1 ) need not be optimal. 

Recall that i 1 T 1 = e 3 — v 2 . Finally, we expand the direct path costs similarly and compute: 



2 (v 2 - v i )+2 2 (e 3 - v 2 ) - Io(v 2 - v 1 ) + X 2 (e 3 - i 1 ) 



> 1 T 1 
~ 2 

1 

_ 2 



~-r 2 yl-e ) - {-r 2 y\ - 6 f 



+ -T 



2^2 



Ml! , Ml 

rpl "T ji 2 



+ 2r 2 vf 



— 77 — On 

rp2 U 



Mo) 2 



Since r 2 > and w 2 can assumed to be positive, the last term is strictly positive, establishing the result. (If 



v\ = then the theorem holds trivially by convexity.) 



□ 



In general, we apply this result under the condition that ri,r 2 > 0. By symmetry it is easily seen that 
the result applied to rotational variations of the paths involved in the result. 

Lemma 15. Consider an RSVP with v, v <G int(M.+), where v = (0,v\,v 2 ) and v = (vi,0, v 2 ). If r\ > r 2 , 
v 2 > vi, then 

ii(v)>i 2 {v). 
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Proof. Let (x* (t) , y* (t) , z* (t)) be an optimal triple corresponding to T\(v), and let T* be the corresponding 
optimal time. It is clear that z*(t) and y*(t) are constant functions due to the convexity property of Lemma 
Thus, we set z* := z*(t) = (0,4,4)' and V* = = (2/1, >°)'- Therefore, x*(t) = z* - Ry* and z* 
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satisfies zJT* = V\, and z|T* = v^- It is clear that z\ < z 2 since «i < v^- 

Setting z(t) = z = (2*,0, V(t) = y= (0,y|,0)' and x(t) = I- Ry, we note that (x(t),y(t), z(t)) is 
a feasible one-piece triple from the origin to v' where z(-) € and T* is the corresponding time for this 
path to reach Then 

Mv) = \\\z*{t) - Ry*(t) - 9\\ 2 T* =\\\z*- Ry* - 9\\ 2 T* . 
On the other hand, since (x(t),y(t), z(t)) is feasible, 

Mv) < H mMt) (v) = l -\\z-Ry-8\\ 2 T*. 

So 

Xx(v) - Mv) > l -(\\z* - Ry* - 6\\ 2 - \\z - Ry - 0\\ 2 )T* = 2(n - r 2 )( 7o - 7 i)(4 - 4MT* > 0. 

The last inequality follows from our assumptions and because 70 > 71, due to Lemma |4| □ 

Our study of optimal path characterizations rests crucially on comparing the paths depicted in Figure [3] 
For paths with both a finite number of segments and an infinite number of segments, we wish to establish 
that the blue path is "cheaper" than the red path. In most instances it is difficult to establish this as 
a general property, so we provide sufficient conditions for this blue-path-red-path condition to hold. For 
specific numerical instances of an RSVP, these conditions can verified easily. Furthermore, a combination of 
numerical and analytical arguments can be used to show that the conditions hold in general on Rf, defined 
below. 

Let Rf = {(ri,r 2 ) € H+ | n > r 2 ,-l <r 1 +r 2 < 2}. 

Condition 1. For an RSVP with reflection matrix R, (ri,f2) S Rf and 

(1 + r|)(l + r\ - r 2 - r x r 2 ) 2 > 2(Yir 2 ) 2 (l + r\ + r\ - r x - r 2 - r x r 2 ). 

Condition 2. For an RSVP with reflection matrix R, (ri,r 2 ) € Rf and at least one of the two inequalities 
below hold for each (vi, v 2 , W3 

[(1 + n + r 2 )T* + r lV2 + r 2 (v 3 - l)] 2 - (1 + rj + r 2 )(2T*i; 1 - vf) > 0. 



(1 + n + r 2 )T* + ^{ n v 2 + r 2 (v 3 - 1)) 



(1 + r{ + r 2 ) ( 9(T*) 2 + 2T* Vl -v\- \\v\ + (v 3 - l) 2 ] ) > 0, 



where 



T = j v 2 + v 2 2 + (v 3 -iy ^ 
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Is the blue path cheaper? 



Figure 3: Red Path - Blue Path Comparison 

These two conditions are required to prove the next two results. Lemma [16] is proved in the Appendix. 

Lemma 16. Given an RSVP with T = I, 9q < 0, and r\ > r-i > 0, define points v — (ui, 0,1)3) and 
v 1 = (0, i»2) V3) with Vi > for i = 1,2, 3. Suppose further that V2 < V3 and vi < 1. Then 

(a) If Condition 1 holds, thenX^iv) > 22(63). 

(b) If Condition 2 holds, then Xq(v' — v) > I\{v' — e^). 

Lemma 17. Given an RSVP with T = I, 6q < 0, and r\,r% > 0, define points v — (wi,0, V3) and v' = 
(0, t>2, V3) with Vi > for i = 1, 2, 3. Then if Conditions 1 and 2 hold, the least cost two-piece path from the 
origin to v and from v to v' is not an optimal path to v' . 



Proof. Invoking the scaling property of Lemma 11 we assume without loss of generality that vs = 1. When 



V2 > V3, the two-piece path through v is not optimal due to Theorem 12 

Case 1. Consider then the case where v-i < v 3 and r 2 > r\ and set v — (0, t>i, 1). We claim that the 
optimal two-piece path through v is strictly better than the optimal two-piece path through v. First, Lemma 
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gives 22 (v) > Ti(v). In other words, the first segment of the path through v has a lower (or equal) cost 
than the first segment through v'. 

To compare the second segments, note that 

2o(i/ — v) = || #|| \\v' — v\\ — (9, v — v) and 
To{v' -v) = \\6\\\\v' -v\\ - (8,v' -v). 
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Furthermore, we have (9, v' — v) = (9, v' — v), and 



\\v' - v\\ = ^Jv 2 + v 2 + (v 3 -l) 2 > ^{v 2 - t'i) 2 + (v 3 - l) 2 = \\v' - v\\. 

Thus, Xq(v' — v) > Io(v' — v) and the result is established for this case. 

Case 2. Suppose next that v 2 < v 3 and v\ > 1 (with no restriction on r\ and r 2 ). Let v = (0,1, Vi) 
and consider the two-piece path to v' via v. Again, we show that the optimal two-piece path through v is 
strictly better than the optimal two-piece path through v. By rotational symmetry T 2 (v) — T\(v). On the 
other hand (9, v' — v) = (9, v' — v), and 

- v|| = ^Jvf + v% + (v 3 -l) 2 > y/(v 2 -l) 2 + (v 3 -v 1 ) 2 = \\v' - v\\. 

As in Case 1, this implies 

2 {v' —v)> 1q(v - v), 

which establishes the result for this case. 

Case 3. The remaining case is when v% < Va, r% < n, and v\ < 1. Once again we find an alternate two- 
piece path to v' which has a lower cost. In this case, consider the two-piece path to v' via e 3 . If Conditions 
1 and 2 hold, then Lemma 16 indicates that > ^2(^3) and io( v ' ~ v ) ^ i\W ~ e 3 ). This establishes 



the result for this case. □ 

Note that Conditions 1 and 2 are only needed to establish the third case. It may be possible to replace 
these conditions by simpler expressions for special cases. 

Next, we are now able to establish the result that there always exist gradual optimal paths to points on 
the boundary of the octant. It is important to note that the class of gradual paths do not include paths 
which traverse an axis and then cross the interior to a point on a two-dimensional face. 

Theorem 18. Consider an RSVP with V = I, 9q < 0, r\,r2 > 0. Suppose Conditions 1 and 2 hold and 
that there exists an optimal path with a finite number of segments. Then: 

(i) For any point on an axis there exists an optimal path consisting of a single segment; and 

(ii) For any point on a two-dimensional face there exists an optimal gradual path, consisting of two seg- 
ments. 

Proof. To prove the result we need to eliminate a large number of path types. In order to categorize these 
types, note that each type can be classified according to the endpoints of the linear segments. The endpoint 
of each piece can be on the interior of a two-dimensional face (F), on the interior of an axis (A), or the 
origin (O). In all the arguments below, we consider a path with a finite number of pieces, and thus a finite 
number of endpoints, which starts at a point v in the octant and terminates at the origin. Specifically, we 
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label the endpoints in "reverse order." Note that an endpoint cannot be in the interior of the octant due 



to the convexity property in Lemma 11 Furthermore, note that the last point is always of type O and of 
course, this is the only position at which this type occurs. 

Next, consider an endpoint of type F. There are three possibilities for the previous endpoint: 

• The endpoint is on an axis A (either one of the two axes adjoining this face, or the remaining axis) 

• The endpoint is on the same face SF. 

• The point is on a difference face DF. 

Similarly, for an endpoint of type A, there are two possibilities for the previous endpoint: 

• The point is on an axis A. 

• The point is on the same face SF (i.e., a face adjoining the axis). 

Notice that for a point of type A the previous point cannot be on the face not adjoining the axis as a 
consequence of the Bad Faces Theorem. With this notation, we can categorize a piecewise linear path by a 
finite scries whose elements arc in the set {SF, DF, A, O}. 

For a series corresponding to a finite-piece optimal path, we infer the following rules: 

1. By the convexity property of optimal paths, none of the following pairs can appear in the series: 
(DF, DF), (SF,SF), (SF,A). 

2. If A appears somewhere in the series, then the end of the series cannot be (A,0), due to the scaling 
and symmetry properties of optimal paths. The only exception is, of course a series which is simply 
(A,0). 

3. The series cannot end with (SF, O) by convexity. 



4. The series cannot end with (DF, O) due to Lemma 17 



Note that rules 3 and 4 imply that the series must end with (A, O) . 

We now establish part (i). Consider a path with the terminal point on say axis Fi t 2 and the first segment 
of the path emanating from the origin. If this first segment traverses an axis, then by scaling and symmetry, 
part (i) immediately holds for any terminal point on an axis. By convexity, the first segment cannot be in 
the interior of the octant. So, the first segment must be embedded in a two-dimensional face. Now, the 
second segment cannot be embedded in this same face due to convexity. So, it must cross the interior and 
terminate either in a different face, or on the opposing axis. The first case is ruled out by Lemma |17| The 
second case is not possible by the Bad Faces theorem. Hence, the first, and only segment, must be embedded 
in an axis. 
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In consideration of Rules 1 through 4 above, to prove part (ii) we must exclude two remaining cases: 
(F, SDFi , DF, A, O) and (F, DSF t , DF, A, O), i = 0, 1, 2, . . . . Here SDF, is a subsequence of (SF, DF) that 
repeats i times, and DSFi is a subsequence of (DF, SF) that repeats i times. Consider the case (F, DF, A, O) 
first. Without loss of generality let the terminal point be on face F\, denote it v 1 — (0, v\,v\), and assume 
that < v\ < v\. The endpoint before v 1 has to be a point v 2 on F^. So v 2 = (v\, 0, v 2 ) which is the DF in 
the series. We must have v\ < v 2 for the path to be optimal, by the assumed type of path and the Bad Faces 
Theorem. The next endpoint v 3 cannot be on axis F 1 3 again by the Bad Faces Theorem. Furthermore, it 

3- 



cannot be on axis F\ 2 due to the arguments in the proof of Lemma 17 Hence, v 3 must be in F2 



Now, if the path just described is optimal, this implies that the optimal path to v 2 is from the origin to 
v 3 then to v 2 . Then by the scaling property, the optimal path to an arbitrary point (zti,0,M3) on face F 2 is 
of the form (F,A,0) if u\ < U3. By symmetry, the optimal path to an arbitrary point (0,^2,^3) on F± is 
also of the form (F,A,0) if U2 < u%. However, v 1 is indeed of this form, which means we can replace the 
proposed optimal path of the form (F, DF, A, O) by a gradual path of the form (F, A, O). 

Next consider the case (F, SF, DF, A, O). If the terminal point v is in say Fx then so is the endpoint v 1 
immediately preceding this point. This implies that the optimal path to v 1 is of the form (F, DF, A, O). As 
argued above, we can eliminate this form. All the remaining cases can be eliminate by analogous arguments 
that reduce the end of the series to the (F, DF, A, O) case. We conclude that any optimal path to a point 
on the interior of a two-dimensional face can be reduced to the gradual path forms (F, A, O) or (F, O). This 
establishes part (ii). □ 

Finally, we present the main result for optimal paths with a finite number of segments. 



Theorem 19. Suppose the conditions in Theorem 18 hold for an RSVP. For any point mR^_, if there exists 



an optimal path with a finite number of pieces then there exists a gradual optimal path. 

Proof. Theorem [18] establishes the result for points on the boundary of the octant. By convexity, the last 
segment of an optimal path to an interior point must have an endpoint on the boundary. The result then 



follows directly from Theorem 18 □ 



9 Results for Exotic Paths 

This entire section is devoted to arguing that certain "exotic spirals" cannot be optimal. Depicted in the 
left-hand side of Figure [4] is what we call a classic spiral, a path type which has appeared in other contexts 



in the literature on fluid models. In Section 10 we show that such a path can indeed be the optimal solution 
to the variational problem we consider in the paper. For now, however, we wish to show that other types of 
paths, exotic spirals, cannot be better than a classic spiral. One important type of exotic spiral appears in 
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Figure 4: Spiral paths 

the right-hand side of Figure |4j Eliminating this type of path from consideration is the focus of much of the 
next two results. 

Lemma 20. For any two-dimensional face, define the bisecting ray to be the ray which forms an angle of 
7r/4 radians with the adjacent axes and whose endpoint is the origin. Consider an optimal path which the 
following characteristics: it contains a line segment that intersects the bisecting ray in F\ and it contains 
another line segment that intersects the bisecting ray in F2. Then there exists an optimal path with the 
following characteristics: 

(a) The path has two segments (as defined above) which form the same angle with the bisecting rays (i.e., 
if the segments are rotated to lie in the same face, then they must be parallel). 

(b) The path contains another segment in F3 which intersects the bisecting ray at the same angle. 

(c) The path contains an infinite number of segments. 



Proof. Let v°v 1 6 F\ and v 2 v 3 € F% be the segments which intersect the respective bisecting rays and 
suppose the points are traversed by the optimal path in the order tr, v , v 2 , and v . 

We prove part (a) by contradiction, assuming the segments do not form the same angles with the bisecting 
rays. Consider then the portion of the path from the origin to v . This portion can be rotated and scaled 



to create an optimal path which passes through the point, call it w, where v 2 v 3 intersects the bisecting ray 
in F%- Thus, we can create a path from the origin through w to V3 which is optimal yet has a "kink" at w. 



However, this path cannot be optimal due to reflected convexity (Lemma 23). This establishes (a). 
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It is clear by rotation, scaling, and merging that one can form an optimal path to V3 which intersects 
the bisecting ray in F3. Repeating the process results in the formation of an optimal path with an infinite 
number of such segments. This establishes (b) and (c). □ 

Theorem 21 (Elmination of Exotic Spirals). Consider an RSVP with T = I, 9o < 0, r\,r% > and suppose 
Conditions 1 and 2 hold. For any optimal path to a point on the axis with a countably infinite number of 
segments, there exists another path, with lower or equal cost, which is of the form of the classic spiral (i.e., 
of the form (A, A, A, ... ,0)). 

Proof. We begin with a general principle that holds for paths with an infinite number of segments. Consider 
an optimal path characterization which begins with an A and contains another A at position n, elsewhere 
in the sequence. Then there exists an optimal path whose entire characterization must be identical to the 
(original) characterization starting at position n. This principle follows directly by scaling, rotation, and 
merging and it can be thought of as enforcing a "self-similarity" property of optimal paths. As an example, 
consider an optimal path of the form (A, AS, A, A, A, A, ... , O) where AS is an arbitrary subsequence. The 
principle implies that such a path can be replaced by a classic spiral of the form (A, A, A, ... ,0). 

Now, consider the terminal point of an optimal path, which by assumption lies on an axis and which by 
our convention is represented by the first A in the sequence characterizing this path. If the next endpoint 
lies on an axis, then the path is a classic spiral (or can be replaced by one), based on the principle above. 
So suppose this is not the case. The second endpoint cannot be on a different face (DF) by the Bad Faces 



Theorem (Theorem 12). Thus the only remaining possibility is that the second endpoint is characterized as 
SF, that is, it lies on the interior of one of the two adjoining faces. 

Next, in any place in the sequence only a DF can follow SF by convexity and the Bad Faces Theorem. 
After a DF, either an SF or an A may follow (DF cannot follow, again by convexity). Finally, between 
any two appearances of an A in the sequence, the SF and DF sequences can be assumed to be the same, 
again invoking the self-similarity principle above. Putting all of these observations together, we conclude 
that apart from the classic spiral case, there are only two other general categories of paths with an infinite 
number of segments: 

(i) (A, SDF, O), where (SDF) is an infinite subsequence of (SF, DF). 

(ii) (A, SDF h A, SDF h A, SDF h ...,0),i = 1,2,3,.... 

We now proceed to eliminate these two types of paths. 

Part (i). We consider first a path of type (A, SDF,0). Without loss of generality, assume that the 
terminal point (the first A in the sequence) is e3 and the next pivot point is v 1 € int(F2)- Now the farthest 
face from v 1 can be either F± or F3 (since the point is on F2, this cannot be the farthest face). Suppose v 1 
is strictly closer to F\. Then by the Bad Faces Theorem, the next endpoint must be in F\. However, such a 
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path can be eliminated from consideration by Lemma 14 Therefore we assume that v is closer to F3 than 



Fi and the next endpoint in the path is in F3 again by the Bad Faces Theorem. (If v 1 is equidistant to F\ 
and F3, then the segments to Fi and F 3 have the same costs and we choose the segment going to F 3 ). The 
next point, v 3 is also in the interior of F3 due to our assumption on the path type. Using arguments from 
the proof of the Bad Faces Theorem, it can be shown that v 2 must be closer to F2 than F±. Next, if v 3 is 
closer to F2 than F\ then the resulting path is of "switchback" form. Such a path is suboptimal by Lemma 



13 



Thus, v must be closer to F\ than F 2 - Furthermore, by Lemma 20 e^v 1 is rotationally parallel to 



v 2 v 3 . Hence, after e 3 the faces containing the endpoints are in this order: F2, -F3, F$, F±, Fx, F2, F2, By 

the usual arguments using rotation and scaling, all the DF segments are rotationally parallel. So, for this 
path type, the path is an "exotic spiral" as depicted in the right half of Figure [4] 

In the remainder of the proof, we argue that such exotic spirals can be replaced by classic spirals. So, 



notice that that if 63V 1 is extended to axis ^2,3 , it intersects it at a point we denote by u = (t, 0, 0). Similarly, 
when v 2 v 3 is extended to axis ^2,3) it intersects it at a point w = (fc,0,0). Now, by symmetry, the "shrink 
rate" must be the same on each face. In other words suppose the next two endpoints in Fi are v and v 5 . 



Then if u 4 u 5 is extended to axis Fi.,3, it intersects the axis at the point (0,fc 2 ,0). Furthermore, in order for 
the path to spiral in to the origin, we must have < k < 1. This leads to four cases related to the values of 
t and k. We explore each of these in turn in order to eliminate the corresponding exotic spirals. 

1. < fc = t < 1. In this case, w — u. As the path is assumed to be optimal it is clear that 

loiv 1 -v 2 ) <T 2 {v 1 -u)+l 3 (u- v 2 ). (18) 



Consider any point v 2 £ F3 between v 2 and v 3 and a point v 1 £ F2 between and v 1 such that v 1; u 2 



is parallel to v 1 v 2 . It is possible to choose such points since w = u. Then, using scaling, translation, 



and ( 18 ) we have 

l Q {v l - v 2 ) < ^(v 1 - v 1 ) +l 3 (v 2 - i^+Mv 1 - v 2 ). 

Thus, the segment of the original path from e 3 to v 3 can replaced by the segments through ea, v , v 2 , 
and v 3 . Now, in particular, we can chose V2 = Vs. However, the new path created in this way is strictly 
suboptimal by convexity. Specifically, it can replaced by a lower cost path which goes directly from v 1 
to f 4 . However, we have already argued that a path following this pattern can be eliminated, which 
establishes our argument for this case. 

2. < t < fc < 1. In order to simplify the algebraic manipulations below, define the following costs: 
a = 12(^3 — u), b = I^v 1 — v 2 ), sia = ^(v 1 — u), S2 = is{w — v 2 ), s 2 ka = Is{w — to), where to is the 



extension of v 2 v 3 to axis F\ 3 and s± + s 2 < 1. With this notation the total cost of this exotic spiral 
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path can be written as: 

1 

I(v) = a(l - si) + ^(6 + fca(l - si - s 2 ))fc* = ^— ^- (6 + a(l - si - s 2 )) + as 2 . 
*=o 

The total cost of the classic spiral path with terminal segment uv is: 



»=o 



If the exotic spiral path is strictly optimal then we must have 

(b + a(l — si — s 2 )) + as 2 < 



l-k y v i-r 

which implies 

h / 1 I \ 

(19) 



s<( 1 -*)(r^-i^( 1 -—»)-»)- 



Now consider a new point v 2 on the extension of uv 2 . When v 2 is close enough to v 2 , there exists 
a point v 1 on w 1 e3 such that v 2 v 1 is parallel to v 2 v 1 . The optimal direct costs on the latter two 
segments are proportional and thus we can define a quantity x by the following relation: ^(w 1 — v 2 ) = 
Xoiv 1 — v 2 )(l + x), where x > 0. These segments also define an exotic spiral, similar to the original one 
under consideration, which passes through e^, v 1 , and v 2 . Define f(x) to be the cost of such a path, 
as follows: 

oo 

f(x) - a(l-s 1 (l + x)) + J2( b ^ + x) + ka(l-(s 1 +s 2 )(l + x)))q l 

*=o 

= os 2 (l + x) + Y—q^ 1 + x ) + a ( l - ( s i + s 2)(l + a;))), 

where q = k — (t — k)x. Notice that x parameterizes a set of new exotic spiral paths, with x = 
corresponding to the original path. If there exists an x > such that the corresponding path has 
a strictly lower cost than the original, then the original path cannot be optimal. To show this, it is 
sufficient to show that the derivative of the cost with respect to x is strictly negative at zero. The 
derivative with respect to a; is: 

f , M _ , [b - ajsi + s 2 )](l ~q)-(t- k)(b(l + x)+ <L - (si + S2 )(l + x))) 
{ ' ~ 2 + (1 - q) 2 

Some algebra yields that /'(0) < is equivalent to: 

[b-a(s 1 + s 2 )}(l-k)-(b + a(l-s 1 -s 2 ))(t-k) 
aS2 + (l^W 2 < °' 

which can be rewritten as 

h - < - j-^(*2(l - k) 2 - ( Sl + s 2 )(l - k) - (1 - si - s 2 )(t - fc)). (20) 
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Inequality (20 1 can be established, via (19), if the following holds: 

(«a(l ~ k) 2 - (si + s 2 )(l - k) - (1 - ax - a 2 )(t - k)) > (1 - k)(-^— - — ^-(1 - «i - s 2 ) - s 2 ). 

1 — I 1 — t 1 — K 

This is equivalent to 

s 2 (fc- l)(fe — t) < 0. (21) 
Recall that k < 1 since the path spirals to the origin, t < k by assumption, and s 2 > also by 



assumption. Thus this final inequality is indeed valid, which establishes (20), implying /'(0) < 0. This 
contradicts the assumed optimality of the exotic spiral. 

3. < k < 1 < t. The arguments for this case are analogous. In particular, the expression for /'(0) is the 
same: 

[6 - a( Sl + s 2 )](l - k) - (b + o(l - Sl - s 2 ))(t - k) 

/(0) = aS2 + a^P 

, [6-a(«i + « 2 )](l-*)-(6 + o(l-si-s 2 ))(l-*) 

- aS2 + a^fc? 

a 

= aS2 -l^k- 

Since s 2 and k are positive, and strictly smaller than 1, the last expression is strictly negative. Thus, 
as in the previous case we conclude that the original exotic spiral cannot be optimal. 



4. < k < t < 1. We use the same construction as in Case 2, note however that (21 ) no longer holds, so 
a modified argument is necessary. So, again suppose that the given exotic spiral path is optimal. As 
argued previously, we must have f'(0) = for the path to be optimal, i.e., 

[6 - a (si + s 2 )](l - fc) - (& + a(l - si - s 2 ))(t - fc) n 

aS2 + a^fc? = 

Rewriting this yields 

b = r ^(-s 2 (l-fc) 2 + ( S i + S2 )(l-fc) + (l- Sl - S2 )(t-fc)) 
asi(l — k) a f M M n 1 . 2 



1-i 1-t 
asi(l — fc) 
1-t 



(s 2 (l - fc) - s 2 (l - k,y + (1 - si - s 2 )(t - k)) (22) 



To prove that the exotic spiral in this case is not optimal, we again propose a path with a strictly lower 
cost. To this end, let v 2 = («i,i>2,0) and define the point v 2 = (vf,0, 0). Condition 2 implies that 

Therefore, if we show that /(w 2 ) < I(y 2 ) then there exists an alternate path with strictly lower cost. 
Using symmetry and scaling, and the fact that the original path to e% was presumed optimal, we have 
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v\ = fc(l — s 2 ) and 

1 



I(v z ) = fc(l - s 2 )I(v) = k(l - s 2 ) 



Furthermore, I(v 2 ) = jzrk(b + a(l — Si — s 2 )). So 



| (b + a(l - si - s 2 )) + as 2 



I(v 2 ) I(v 2 ) = ks 2 + b - r ^ + r A_a(l - S2 ) > 0, 



which is the desired result. Note that the second term is strictly positive due to (22). So, again, this 
contradicts the optimality of the exotic spiral. 

Part (ii). We now turn our attention to the other general type of exotic spiral, one with the character- 
ization (A, SDFi, A, SDFi, A, SDFi, ...,0), i € {1,2,3,.. .} where each SDFi is a sequence with SF/DF 
segments repeated i times. By symmetry and scaling arguments, we can assume that each of these sequences 
is identical. Suppose that the end point v° is on axis i*\ 2 and the SF is from v 1 on F\. By Lemma 
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v 2 is on a different face than v 1 , v 1 must be closer to than to F\^ 2 and v 2 must be in F3. The next point, 
v 3 , is either on axis F23 or in F3. However, it must be closer to -F2.3 than to i*i 3. Then based on Lemma 
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there exists point u on axis Fx 3 for which an optimal path to u contains v 3 v 2 . By scaling and symmetry, 
this path can be assumed to be of the (A, SDFi, A, SDFi, A, SDFi, . . . , O) form posited for the path to i>o- 
Now, by assumption, there are i SF/DF segments between v° and the next point on the axis and of course 
this path passes through v 2 . Considering the optimal path to u, since it also passes through v 2 en route to 
u, the portion of the this path to v 2 can be replaced by the optimal path to Vq, up to v 2 . This patching 
process forms another optimal path to u. Since there are i SF/DF segments between v° and the next axis 
point, there are then i — 1 SF/DF segments between u and this same point. Therefore we have constructed 
an optimal path to a point on the axis (u) which is of the form (A, SDFi-i, A, SDFi^i, A, SDFi_i, . . . , O). 
Hence, there must exist an optimal path of the same form to v . Repeating this patching process results in 
the construction of an optimal path to (and thus any point on any axis) of the form (A, A, A, ... ,0). 

It remains to be argued that there must be a finite number of segments between in an optimal path 
between any two axis points. We only give an outline here. Consider a path with an infinite number of 
segments which converge to a point vq on i*i 2- There must be an infinite number endpoints of such segments 



in an e-ball around vq- Furthermore, by Lemma 13 there must exist an infinite subsequence of endpoints 



for which the other terminal point of the segment is in F3. The cost of all such segments can be uniformly 
bounded away from zero (using the infimum of the cost from the e-ball to F 3 , which is strictly positive). 
This implies however, that the total cost of any such path is infinite. Hence, the path cannot be optimal. □ 
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10 An Optimal Spiral Path 

In Example 2 of Section 6 of El Kharoubi et al. [5], it is shown that a spiral path has a lower cost than a 
two-piece gradual path, for the corresponding RSVP. Here we give a related example, and using the results 
of previous sections, show that a spiral path is indeed optimal. 
Let 6 = (-1, -1, -1)', r = I and 

/ 1 | \ 
R= I 1 ■ 

V 1 W 

To establish that a spiral path is optimal, we need to undertake three steps. First, we check that Conditions 
1 and 2 of Section [8] hold. If so, then we know that only gradual paths or spiral paths are optimal. Second, 
using results from [5] we check the reflectivity characteristics of optimal paths traversing an axis. Third, to 
travel to a point, say e^, on an axis, we verify that it is less costly to traverse one of the other axes and then 
cross a two-dimensional face. If this is the case, then one can construct a spiral path to e 3 that is cheaper 
than the gradual path to e% (which simply travels along the axis). 

The next proposition simplifies the process of checking Conditions 1 and 2 and may be useful in producing 
other examples. 

Proposition 22. For a RSVP with r 2 = and < r\ < 2 Conditions 1 and 2 hold. 

Proof. Clearly (ri,r 2 ) € Rf under the assumptions given. Recall that Condition 1 is given by 

(1 + r§)(l + rj - r 2 - nr 2 ) 2 > 2(nr 2 ) 2 (l + rj + r\ - n - r 2 - nr 2 ). (23) 

For r 2 — the condition reduces to (1 + r 2 ) 2 > which clearly holds for all real r\. 

Next, to check Condition 2, it is sufficient to establish that for all v — (vi, 0, 1) G F 2 and v' = (0, v 2 , V3) € 
Fi the following expression is non-negative: 

[(1 + n + r 2 )T* + r lV2 + r 2 (v 3 - l)] 2 - (1 + r\ + r 2 2 ){2T*v 1 - v\). (24) 

By basic algebra, for any real T* and v\ we have 2T*v\ — vf < (T*) 2 . Using this, plus r\ > and r 2 = 



yields the following (non-strict) lower bound on (24): 



[(1 + n)T*} 2 - (1 + r\){T*f = (T*) 2 [(l + n) a - (1 + r\)\ > 0. 
Therefore, Condition 2 also holds. □ 

So, the proposition provides verification of Conditions 1 and 2 for the example in this section. Next, we 
use results from [8] to check reflectivity of the axes. In particular we use equations (24) and (25), Remark 
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2, and Proposition 1 from that paper. Let R\ = (1,3/2,0)', R 2 = (0,1,3/2)', and Ri j2 = (2?i,i?2)- Define 
Ai.2 = I — RxpBxp and B lt2 = {R\ 2-^i,2) _1 -Ri 2- Some algebra shows that 

, 7 1,a ^ . .Bi,2e3 - B lj2 e « (0.0526, 1.5526) > 0. 

So optimal one-piece reflected paths confined to an axis use both corresponding reflection vectors. 
Finally, we check the spiral condition for the point e 3 . In particular either 

24,2 (e 3 ) >!{\ 3 },2( e 3) or Xi,2(e 3 ) >I { 2 i !3}i i(e 3 ) 
must hold. We verify the first inequality. For the parameters of our example we have 

Zi,a(e 3 ) = \\A h2 0\\\\A h2 e 3 \\ - (A h2 6,A h2 e 3 ) « 0.4211. 
Next, let u = (0.5,0,0). Then 

^{2, 3 },2( e 3) < 2^2,3 («) +2: 2 (e 3 - it) ps 0.3317 < 2r 1)2 (e 3 ). 

Therefore, for the given RSVP, the optimal path to any point on the boundary of the octant is a classic 
spiral optimal path. 

We can more precisely characterize this optimal spiral path. In particular, the last piece connects the 
points kei and e 3 , where fe, < k < 1, is the shrink factor. The optimal value of k can be calculated by 
defining the corresponding spiral cost as a function of k: 

f(k)= i ^: k k ei) . 

Applying the data for this problem and setting f'(k) = results in the quadratic root-finding problem 
1228123A: 2 - 3690960A: + 1626300 = 0. The appropriate root is k* ps 0.5363. From this we can calculate the 
cost of the optimal spiral as 

2^2 (e 3 — fc*ei ) 

— ^ re 0.2384. 

1 - k* 

11 Conclusions 

As mentioned in the introduction, this paper only provides a piece of the puzzle of the variational problem 
related to large deviations for SRBM in the orthant. We have only addressed problems with symmetric 
data, and even then some of the results have further restrictions on the parameters. Although we do not 
provide an analytical proof, in Liang |12j convincing numerical evidence indicates that Conditions 1 and 
2 hold whenever the SRBM is stable and ri,r2 > 0. Other results require that the covariance matrix T 
is the identity. This is one condition that seems difficult to remove mathematically, despite the fact that 
intuition indicates that the covariance matrix should not affect the types of paths which are optimal for a 
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given SRBM. We believe that the restrictions on R and T can be removed with some effort, but as noted 
earlier, it is already know that our results cannot be generalized to arbitrary (stable) SRBM data. The 
example in [7 a which was discussed in the introduction is particularly troubling because the reflection matrix 
in that example is partially rotationally symmetric (two of the three reflection vectors exhibit rotationally 
symmetry) . 

In addition, even once one has a handle on the types of paths which are optimal, computation and 
comparison of these path costs appear still requires considerable effort. This indicates that fully solving 
large deviations problems in high dimensions is likely to remain a challenge and future work in dimensions 
four and higher will require examination of specialized cases and increasing mathematical creativity. 
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12 Appendix 

Proof of Lemma^ It can be checked that if T _1 has the form of the lemma, then for a rotationally symmetric 
matrix T the nine equations in FF _1 = I are consistent. Since matrix inverses are unique, it follows 
immediately that F _1 can be written as stated in the lemma. Next, in order for rr _1 = I to hold we must 
have 

7o + 2p7! = 1 and 
7i+p(7o+7i) = 0. 



Solving these equations yields 



which implies 



7i 1 -7o / otn 

P= ; = ~z 1 25 ) 

7o + 7i 2 7i 



2 7l 2 = ( 7o -l)( 7 o+7i). (26) 
Note that 70 = is not possible. To prove 70 > 71 we examine four cases. 
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1. If 70 > and 71 < then the result follows immediately. 

2. Suppose 70 > 0,71 > and 70 < 7i- Then we have 

(70 - l)(7o + 7i) < (71 - l)(7i + 7i) < 2 7i ■ 



This contradicts (26) 



3. Suppose 70 < 0, 71 < and 70 < 7i- Then we have 

(7o - l)(7o + 7i) > (7i ~ !)(7i + 7i) > 2 7i • 



This again contradicts ( 26 1 



4. Suppose finally that 70 < 0,71 > and 70 < 71. Since p < 1 by definition, (25) implies that 



-27! + 1 < 7o < 0. Solving (26) gives 



7o 



1 - 71 - V^l! +271 + I ^ 1 - 71 - 371 



-2 7 i + 1/2, 



which contradicts —271 + 1 < 70. In solving (26) we take the smaller root, since we must have 70 < 0. 

Thus, by contradiction, we have established 70 > 7i- □ 

Proof of Lemma^ It can be checked that if has the form of the lemma, then for a rotationally symmetric 
matrix R the nine equations in RRr 1 — I are consistent. Since matrix inverses are unique, it follows 
immediately that can be written as stated in the lemma. Also, from = /, we have 

a + r 2 c + nb = 1 
b + + T\C = 
c + r2& + ria = 0. 



Summing these equations gives the claimed equality. 

Proof of Theorem \10[ Lemma [6] states that R being completely-5 is equivalent to 

n + r 2 + i> 0. 



By definition, the conditions 



1 + rf 



1 - rir 2 > 
! - 3nr 2 > 



□ 

(27) 

(28) 
(29) 



are necessary and sufficient for R to be a P-matrix. We prove that ( 27 ) is equivalent to ( 28 ) and ( 29 ) by 
partitioning the possible values of r\ and r 2 . 
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1. Suppose that either r% = or r 2 = 0. We prove the case in which r 2 = 0, the other case is analogous. 
When r 2 = |27} reduces to n + 1 > 0, (28 1 is trivially satisfied, and ((29l reduces to r\ + 1 > 0. The 



first inequality is equivalent to the last, establishing the result for this case. 



2. Suppose r 1 ,r 2 > 0. It is easy to see that (27 1 always holds in this case. Furthermore, r x + r 2 < 2 
implies 



ri + r 2 



< 1. 



The arithmetic-geometric mean (AGM) inequality gives 



nr 2 < 



ri + r 2 



(30) 



and thus (28 1 always holds. Invoking the AGM inequality again yields 



rir 2 = \l r\r\ < 



Recall that equality in the AGM inequality holds iff the three terms are equal. Equality of the terms 



in this case implies r± = r-2 = 1 which is not possible due to r\ + r 2 < 2. Therefore (29 1 automatically 
holds. 



3. Suppose ri,r 2 < 0. Note that 



1 + r\ + r\ - 3rir 2 = (n + r 2 + 1)(1 + r\ + r\ - r\ - r 2 - rir 2 ) 



(31) 



and 



1 + r? + r\ - n - r 2 - rir 2 > 0, 



when ri,r 2 < 0. Therefore, in this case, (27) and (29) are equivalent and we need only show that (27) 



implies (28). Note then that n + r 2 + 1 > 0, implies that n > — 1 and r 2 > — 1. Given that both n 



and r 2 are also negative this yields r\r 2 < 1. 



4. Suppose ri > and r 2 < (the case ri < 0, r 2 > is analogous). It is obvious that (28) always holds 



in this case and so we need only show that (27) and (29) are equivalent. Consider again the last term 



in (31) 



1 + r\ + r\ - ri - r 2 - r x r 2 = (r 1 - r 2 ) 2 - (n - 1)(1 - r 2 ). 

Note that when t\ < 1 this term is positive as can be seen from the right-hand side above. When 
7"i > 1 we have r% — r 2 > r± — 1 > and r± — r 2 > 1 — r 2 > and again the right-sand side above is 



clearly positive. This fact implies that (27) and (29) are equivalent, as argued in Case 3 



□ 
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Lemma 23 (Reflected Convexity). Consider a section of a feasible triple (x, y, z) in which the path z consists 



of segments v 1 v 2 and v 2 v 3 , with v ,V ,v G Fj where j € {1,2,3}. Suppose that v 1 v 2 is a reflected segment 



and v 2 v 3 is direct. Then there exists a linear reflected path from vi to v 3 whose cost is no greater than the 
original path. 

Proof. Without loss of generality, assume that j = 1. Let (x 1 (t) , y 1 (t) , z 1 (t)) be the triple corresponding to 
the segment v x v 2 with T = T . Similarly, let (x 2 (t) , y 2 (t) , z 2 (t)) be the triple corresponding to v 2 v 3 with 
T = T 2 . Note that (x 1 (t) , y 1 (t) , z 1 (t)) = (a: 1 , y 1 , i x )t and (x 2 {t),y 2 {t), z 2 {t)) = (x 2 ,y 2 ,z 2 )t. Further denote 
x 1 = (x\,x\, x\)' . Similar notation is used for the other variables. Notice that y 1 = (y{,0, 0)' and y 2 = 0. 



By our assumptions on the segments, we have x 1 + Ry 1 



z , x 



0, x\ < and x\ = 0. By 



translation, we set ^(T 1 ) = v 2 - v 1 and z 2 {T 2 ) = x 2 (T 2 ) = v 3 - v 2 . Also, define points u 2 = v 1 + x 1 ^ 1 ) 
and u 3 = u 2 + x 2 (T 2 ). Notice that these two points are not in the interior of the octant. Based on convexity, 



Mu 3 



'') +l (u 2 -v 1 ) >l (u 3 ~ v 1 ) 



Let x 3 (t) be optimal to 1q(u 3 — v 1 ) with corresponding T — T 3 , where x 3 (t) = x 3 t. It is clear that x 3 T 3 = u 3 ~ 
v 1 = x 1 T 1 +x 2 T 2 . Define y 3 = and y 3 (t) = y 3 t. Also define z 3 (t) = z 3 (t) + Ry 3 (t) . So z 3 (T 3 ) = v 3 -v 1 . 
Thus (x 3 (t), y 3 (t), z 3 (t)) is a feasible triple for lj(v 3 - v 1 ). Therefore, 

Ti^-v 1 ) < I f \\x 3 {t)-d\fdt=i Q {u 3 ~v x ) <2 (u 3 - u 2 ) +2 {u 2 - v 1 ) = l^v 2 - v 1 ) +2 Q (v 3 ~ v 2 ), 
2 Jo 



which establishes the result. 



□ 



Proof of Lemma 16 Part (a). The claim is that if Condition 1 holds, then 22{v) > 22(e%). Without loss of 
generality we set v — (i>i,0, 1). For all non-negative v%, define the function 

G(v x ) :=2 2 (v) = \\Av\\\\Ae\\-(Av,A6), 

where A — I — R 2 B, B — (R 2 R 2 )~ 1 R' 2 , and R 2 = (r 2 , l,fi)'. It can be checked that G(-) is strictly convex 



on (0, 1). Therefore, to prove 2 2 (v) > 2 2 (e^) for v\ > 0, it is enough to show that 



9+GQi) 



Ul= o > 0. Some 



algebra yields 



d+G( Vl ) 



dvi 



1 \\A0\ 



.(A 31 + A 13 )-(AB) 1 



2 ||^e 3 |! 

Note that A 3 i + A 13 < and (A9) 1 < in Rf. So, to prove the non- negativity of the derivative, it is 
sufficient to show that 



{A9)l > 



1 



\\A9\\ 
\Ae 3 \ 



(A 31 +A 13 ) 



(32) 



Next, we have 



\Ae 3 \ 



2(1 + r\+rl-r 1 -r 2 - nTg) | 
l + r 2 ' 



'()• 
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and 



Plugging these equalities into ( 32 ) yields the condition 



(1 + rl)(l + r\ - r 2 - nr 2 ) 2 > 2{r 1 r 2 f{\ + r\ + r\ - n - r 2 - nr 2 ). (33) 



In summary, if (33) holds, then & +^( Vl > \ Vl=0 > which in turn implies i 2 {v) > X 2 (e 3 ). 



Part (b). The claim is that if Condition 2 holds then Iq(v' — v) > I\(v' — e 3 ). Since 1q(v' — v) 
and Ii(v' — e 3 ) are both proportional to 6*o it is enough to verify the case when 6*o = —1- We have that 
Z (y -v) = ±\\x*(t) - 9\\ 2 T* where 



mt \\V'-V\\ (Vl) 2 + (V2) 2 + (V3-1) 2 

= ^er = y 3 

and x*(t) = x*t = t(x*, x^, x^)' for < t < T* . We construct a feasible reflected path in F 1 from e 3 to v' with 
a cost Hy.i(v' — 63) less than or equal to T.q(v' — v). This construction then implies Tq(v' — v) > — 63). 
Below, we construct feasible triples with two different "speeds." Of course, in the search for a lower cost 
path, one could examine all possible feasible speeds to generate a weaker sufficient condition than that which 
appears in the lemma, but the cases we consider turn out to be sufficient for our purposes. 

Denote one feasible triple from to v' by (x(t), y(t), z{t)) and T. Let T — T* , z(t) = zt = t(Q, z 2 ,z^)' 
and y(t) = t(y±,Q, 0)' for some y\ > 0. It is clear that z 2 = x 2 and Z3 = £3. The goal now is to determine if 
there exists a y\ > such that 

Io(v' - v) - Hy, 2 (v' - e 3 ) > 0. 

Plugging in x{t) = z{t) — Ry{t) 1 x*T* = v' — v, and writing y — y\T* we see that the inequality above is 
equivalent to 

(T* - v,) 2 + (T* + v 2 ) 2 + (T* +v 3 - I) 2 - [(T* - y) 2 + (T* + v 2 - r x y) 2 + (T* +v 3 -l- r 2 y) 2 ] > 0. (34) 

Denote the second feasible triple by (x(t), y(t), z(t)) and T, where y(t) — t(yi, 0, 0)'. This time however, set 
T = AT*. As in the first case, we have z\ — 0, z 2 — ^ and £3 = Let y = y\T* . The goal again is to 
determine if there exists a y\ > such that 

l {v' - v) - H v ,z(v' - e 3 ) > 0, 

which is equivalent to 

(T* - v,) 2 + (T* + v 2 ) 2 + (T* + v 3 - l) 2 - 4[(T* - y) 2 + (T* + | - ny) 2 + (T* + - r 2 y) 2 ] > 0. (35) 

In summary, we claim that for given problem data and points v and v' if there exists ay>0ory>0 



such that (34) or (35), respectively, is satisfied, then the desired feasible path construction can be achieved. 
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Notice that the left-hand sides (LHS) in these two equations are both concave, quadratic functions of y and 
y, respectively. So to prove the desired inequalities, it is necessary that at least one of these two functions 
has a non-negative maximum which is achieved at a non-negative value. 



In (34), the maximum value of the LHS is reached when 



_* _ (1 + ri + r 2 )T* + r±V2 + r 2 (v 3 - 1) 
V l + rl + rl 

and this maximum is achieved at a non-negative value when 



[(1 + n + r 2 )T* + nv 2 + r 2 (v 3 - l)] 2 - (1 + r\ + r 2 2 ){2T* Vl - v{) > 0. 



(36) 



When U3 > 1 it is easy to see that y* > 0. Considering now the LHS of (36) we have that 



LHS > (1 + n + r 2 y(T*Y - (1 + rf + r' 2 ){T*Y = (T*) 2 (2rx + 2r 2 + 2r x r 2 ) > 



So, when U3 > 1, (34) is satisfied. The case when V3 < 1 is more complicated and both (34) and (35) need 



to be checked. In this case we have 

Mj < \t* , 1 / , w (l + r 1 + r 2 )(ui+« 2 + (l-«3)) . , , 1W , . . n 

{l+r 1 +r 2 )l +r 1 v 2 +r 2 {v 3 -l) > ^r 1 v 2 +r 2 {v 3 -l) > r 2 (v 1 +v 2 )+r 1 v 2 > 0. 

o 



So it is always true that the maximum y* is non-negative. Similarly for (35) the maximum value of the LHS 
is reached when 

(1 + n + r 2 )T* + 1 [r lV2 + r 2 (v 3 - 1)] 



V 



1 + r 2 + r\ 



Again, through straightforward algebra, it can be established that y* > 0. The maximum is achieved at a 
non- negative value when 

4[(1 + n + r 2 )T* + l(ri«a + ^ 2 (« 3 - l))] 2 - (1 + r 2 + r 2 )(9(T*) 2 + 2T* Vl - v\ - ? [ W 2 + (« 3 - l) 2 ]) > 0. (37) 



Thus, we conclude that when either (36) or (37) holds we can construct the required feasible path and 
i {v'-v) >li(v'-e 3 ). 



□ 
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